o 

(N 

> 
O 



c 

o 
o 



> 

On 



X 



Energy dissipation and switching delay in stress-induced switching of 
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Switching the magnetization of a shape-anisotropic 2-phase multiferroic nanomagnet with voltage-generated 
stress is known to dissipate very little energy (< 1 aJ for a switching time of ~0.5 ns) at K temperature. 
Here, we show by solving the stochastic Landau-Lifshitz-Gilbert equation that switching can be carried out 
with ^100% probability in less than 1 ns while dissipating less than 2 aJ at room temperature. This makes 
nanomagnetic logic and memory systems, predicated on stress-induced magnetic reversal, one of the most 
energy-efficient computing hardware extant. We also study the dependence of energy dissipation, switching 
delay, and the critical stress needed to switch, on the rate at which stress is ramped up or down. 
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I. INTRODUCTION 

Shape-anisotropic multiferroic nanomagnets, consist- 
ing of magnetostrictive layers elastically coupled with 
piezoelectric layers^, have emerged as attractive storage 
and switching elements for non- volatile memory and logic 
systems since they are potentially very energy-efficient. 
Their magnetizations can be switched in less than 1 
nanosecond with energy dissipation less than 1 aJ, when 
no thermal noise is present 3 -^. This has led to multiple 
logic proposals incorporating these systems^—. The mag- 
netization of the magnet has two (mutually anti-parallel) 
stable states along the easy axis that encode the binary 
bits and 1. The magnetization is flipped from one 
stable state to the other by applying a tiny voltage of 
few tens of millivolts across the piezoelectric layer while 
constraining the magnetostrictive layer from expanding 
or contracting along its in-plane hard-axis. The voltage 
generates a strain in the piezoelectric layer, which is then 
transferred to the magnetostrictive layer. That produces 
a uniaxial stress in the magnetostrictive layer along its 
easy-axis and rotates the magnetization towards the in- 
plane hard axis as long as the product of the stress and 
the magnetostrictive coefficient is negative. Large angle 
rotations by this method have been demonstrated in re- 
cent experiments^, although not in nanoscale. 

In this paper, we have studied the switching dynam- 
ics of a single-domain magnetostrictive nanomagnet, sub- 
jected to uniaxial stress, in the presence of thermal fluc- 
tuations. The dynamics is governed by the stochas- 
tic Landau-Lifshitz-Gilbert (LLG) equation*^ that de- 
scribes the time-evolution of the magnetization vec- 
tor's orientation under various torques. There are three 
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torques to consider here: torque due to shape anisotropy, 
torque due to stress, and the torque associated with ran- 
dom thermal fluctuations. With realistic ramp rates 
(rate at which stress on the magnet is ramped up or 
down) a magnet can be switched with ^100% probability 
with a (thermally averaged) switching delay of ^0.5 ns 
and (thermally averaged) energy dissipation ^200 kT at 
room-temperature. This is very promising for "beyond- 
Moore's law" ultra-low-energy computing 11-13 . Our sim- 
ulation results show the following: (1) a fast ramp and a 
sufficiently high stress are required to switch the magnet 
with high probability in the presence of thermal noise, 
(2) the stress needed to switch with a given probability 
increases with decreasing ramp rate, (3) if the ramp rate 
is too slow, then the switching probability may never ap- 
proach 100% no matter how much stress is applied, (4) 
the switching probability increases monotonically with 
stress and saturates at ^100% when the ramp is fast, but 
exhibits a non-monotonic dependence on stress when the 
ramp is slow, and (5) the thermal averages of the switch- 
ing delay and energy dissipation are nearly independent 
of the ramp rate if we always switch with the critical 
stress, which is the minimum value of stress needed to 
switch with non-zero probability in the presence of noise. 



II. MODEL 

A. Magnetization dynamics of a magnetostrictive 
nanomagnet in the presence of thermal noise: Solution of 
the stochastic Landau-Lifshitz-Gilbert equation 

Consider an isolated nanomagnet in the shape of an 
elliptical cylinder whose elliptical cross section lies in the 
y-z plane with its major axis aligned along the z-direction 
and minor axis along the y-direction (Fig.[T]) The dimen- 
sion of the major axis is a, that of the minor axis is b, 
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FIG. 1. A two-phase multiferroic nanomagnet in the shape of 
an elliptical cylinder is stressed with an applied voltage via the 
c?3i coupling in the piezoelectric. Mechanical constraints (not 
shown) prevent expansion and contraction of the multiferroic 
along the in- plane hard axis (y-axis). 



and the thickness is I. The z-axis is the easy axis, the 
y-axis is the in-plane hard axis and the x-axis is the out- 
of-plane hard axis. Since I <C b, the out-of-plane hard 
axis is much harder than the in-plane hard axis. Let 9(t) 
be the polar angle and (j>(t) the azimuthal angle of the 
magnetization vector. 

The total energy of the single-domain, magnetostric- 
tive, polycrystalline nanomagnet, subjected to uniaxial 
stress along the easy axis (major axis of the ellipse) is 
the sum of the uniaxial shape anisotropy energy and the 
uniaxial stress anisotropy energy^. The former is given 
byii E SHA (t) = (^ Q /2)M 2 riN d (t), where M s is the sat- 
uration magnetization and Nd(t) is the demagnetization 
factor expressed asii 

N d (t) = N d - zz cos 2 6{t) + N d - yy sin 2 6(t) sin 2 <j){t) 

+ N d ^ xx sin 2 9{t)cos 2 <j)(t) (1) 

with N d - zz , N d - yy , and N d ^ xx being the components 
of the demagnetization factor along the z-axis, y-axis, 
and rr-axis, respectively-^. These factors depend on the 
dimensions of the magnet (values of a, & and I). We 
choose these dimensions as a = 100 nm, b = 90 nm 
and I — 6 nm, which ensures that the magnet has a 
single domain^. These dimensions also determine the 
shape anisotropy energy barriers. The in-plane barrier 
Eb, which is the difference between the shape anisotropy 
energies when 9 = 90° and 9 = 0°, 180° {<j> = ±90°) de- 
termines the static error probability, which is the proba- 
bility of spontaneous magnetization reversal due to ther- 
mal noise. This probability is exp [— E^/kT]. For the 
dimensions and material chosen, Eb = 44 kT at room 
temperature, so that the static error probability at room 
temperature is e -44 . 

The stress anisotropy energy is given byi 4 - EsTA(t) = 
~(3/2)A s <r(t)ilcos 2 9(t), where (3/2)A s is the magne- 
tostriction coefficient of the nanomagnet and a(t) is the 
stress at an instant of time t. Note that a positive X s c(t) 
product will favor alignment of the magnetization along 
the major axis (z-axis), while a negative X s a(t) product 
will favor alignment along the minor axis (y-axis), be- 
cause that will minimize EsTA(t)- In our convention, a 
compressive stress is negative and tensile stress is posi- 
tive. Therefore, in a material like Terfenol-D that has 



positive A s , a compressive stress will favor alignment 
along the minor axis (in-plane hard axis), and tensile 
along the major axis (easy axis)£. 

At any instant of time t, the total energy of the nano- 
magnet can be expressed as 

E(t) = E{9{t), (f>(t), a{t)) = B{t)sin 2 9{t) + C{t) (2) 

where 

B(t) = B a (t) + B stress (t) (3a) 
B (t) = { m /2)M 2 n[N d - xx co S 2 c^{t) 

+ Nd-yysin 2 ^) - N d ^ zz ] (3b) 

B s tress(t) = (3/2)A,ff(t)n (3C) 

C(t) = {fi /2)M 2 nN d ^ zz - (3/2)A s a(t)Q. (3d) 

The torque acting on the magnetization per unit vol- 
ume due to shape and stress anisotropy is 

T E (i) = -n m (t) x VE(e(t),<t>(t),a(t)) 

= -2B(t)sin9(t)cos9(t) e - B Qe (t)sin9(t) e e (4) 

where B 0e (t) = (/i /2) M 2 n{N d ^ xx - N d ^ vv )sin(2<j){t)). 

The torque due to thermal fluctuations is treated via 
a random magnetic field h(i) and is expressed as 



h(t) = h x (t)e x + h y (t)e y + h z (t)e z 



(5) 



where h x (t), h y (t), and h z (t) are the three compo- 
nents of the random thermal field h(t) in x-, y-, and 
z-direction, respectively in Cartesian coordinates. We as- 
sume the properties of the random field h(t) as described 
in Ref. [10j . Accordingly, the random thermal field can 
be expressed as 



hi(t) 



2akT 



|7|(l + a 2 )M y At 



G (os){t) (i = x,y,z) 



(6) 

where 1/Ai is proportional to the attempt frequency of 
the thermal field. Consequently, At should be the sim- 
ulation time-step used to simulate switching trajectories 
in the presence of random thermal torque. The quantity 
G(n t i\{t) is a Gaussian distribution with zero mean and 
unit variance 1 ^. 

The thermal torque can be written as 

TthO) = My n m (i) x h(t) = Pg{t) e> - P (i) e g (7) 

where 

P e (t) = M v [h x {t) cos9{t) cos(f>{t) + h y (t) cos6(t)sin(j>(t) 
-h z (t) sin9{t)\ (8) 
P (i) = M v [h y (t) cos(j>{t) ~ h x (t) sincj>(t)]. (9) 

The magnetization dynamics under the action of the 
torques Te(£) and Tth(^) is described by the stochastic 
Landau-Lifshitz-Gilbert (LLG) equation as follows. 



dn m (t) ( dn m (t) 
-a n m (t) x 



df 



df 



-^[Te^ + TthM] 
(10) 
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where a is the dimensionless phenomenological Gilbert 
damping constant, 7 = 2/xbMo/^ is the gyromag- 
netic ratio for electrons and is equal to 2.21 x 10 s 
(rad.m).(A.s) -1 , /j,b is the Bohr magneton, and My = 

From the last equation, we get the following coupled 
equations for the dynamics of 9(t) and (j>(t). 



(1 



a 



dOit) | 7 | 



dt 



±LL[B 0e (t)sin6(t) 



2aB(t)sin9(t)cos9(t) + (aP s + P )]. (11) 



(1 + a 



#(t) 



Ji[aS e(i)+2S(i)cos&(«) 

My 



- sin~ l 9(t)(P g - aP$)]. (sine ^ 0.) (12) 

These equations describe the magnetization dynamics, 
namely the temporal evolution of the magnetization vec- 
tor's orientation, in the presence of thermal noise. 



B. Fluctuation of magnetization around the easy axis 
(stable orientation) due to thermal noise 

The torque on the magnetization vector due to shape 
and stress anisotropy vanishes when sinO = [see Equa- 
tion i.e. when the magnetization vector is aligned 
along the easy axis. That is why 9 = 0°, 180° are called 
stagnation points. Only thermal fluctuations can budge 
the magnetization vector from the easy axis. To see this, 
consider the situation when 9 = 180°. We get: 



<j)(t) = tan 



ah y (t) + h x (t) 
ah x (t) 



hy(t) 



(13) 



0'(*) = -|7l 



hl(t)+hl(t) 



y/(h y (t) ~ah x (t)) 2 + ( ah y (t) + /l x (i)) 2 



(14) 

We can see clearly from the above equation that ther- 
mal torque can deflect the magnetization from the easy 
axis since the time rate of change of 9(t) [i.e., 9'(t) is 
non-zero in the presence of the thermal field. Note that 
the initial deflection from the easy axis due to the ther- 
mal torque does not depend on the component of the 
random thermal field along the z-axis, i.e., h z (t), which 
is a consequence of having ±z-axis as the easy axes of 
the nanomagnet. However, once the magnetization di- 
rection is even slightly deflected from the easy axis, all 
three components of the random thermal field along the 
x-, y-, and z-direction would come into play and affect 
the deflection. 



C. Thermal distribution of the initial orientation of the 
magnetization vector 

The thermal distributions of 9 and <j) in the unstressed 
magnet are found by solving the Equations (fTTj) and ([12")) 



while setting B stre ss — 0. This will yield the distri- 
bution of the magnetization vector's initial orientation 
when stress is turned on. The ^-distribution is Boltz- 
mann peaked at 9 = 0° or 180°, while the (^-distribution 
is Gaussian peaked at 4> = ±90° (Ref. [18| ). Since the 
most probable value of 9 is either 0° or 180°, where stress 
is ineffective (stagnation point), there are long tails in 
the switching delay distribution. They are due to the 
fact that when we start out from 9 = 0°, 180°, we have 
to wait a while before thermal kick sets the switching 
in motion. Thus, switching trajectories initiating from a 
stagnation point are very slo w 19 i 20 . 

In order to eliminate the long tails in the switching de- 
lay distribution and thus decrease the mean switching de- 
lay, one can apply a small static bias magnetic field that 
will shift the peak of 9i n iu a i distribution away from the 
easy axis, so that the most probable starting orientation 
will no longer be a stagnation point. This field is applied 
along the out-of-plane hard axis (+x-direction) so that 
the potential energy due to the applied magnetic field be- 
comes E mag (t) = —MvHsin9(t) cos<f>(t), where H is the 
magnitude of magnetic field. The torque generated due 
to this field is Tm (t) = -n m (t)xVE mag (e(t),cj)(t)). The 
presence of this field will modify Equations (JTTJ) and (fT"2"|) 
to 



(1 + a 2 ) = -j^-[B e(t)sin9(t)-2aB(t)sin9(t)cos9(t) 

+ aM v H cos9(t) cos<j>(t) - M V H sin<j)(t) 

+ (aP e + P*)}, (15) 



(1 



d<j)(t) 
dt 



j£-[aB 0e (t) + 2B(t)cos9(t) 



~ sin~ 1 9(t) (M V H cos9(t) cos(f>(t) + aM v H sin<j>(t)) 
- sin- 1 9(t)(P e -aP (t> )]. (sin9^0.) (16) 

The bias field also makes the potential energy profile of 
the magnet asymmetric in 0-space and the energy mini- 
mum will be shifted from <\> m i n = ±90° (the plane of the 
magnet) to 



cos 



H 



M s (N d _ a 



(17) 



However, the profile will remain symmetric in #-space, 
with 9 = 0° and 9 = 180° remaining as the minimum en- 
ergy locations. With the parameters used in this paper, 
a bias magnetic field of flux density 40 mT would make 
(f> m in — 87° . Application of the bias magnetic field will 
also reduce the in-plane shape anisotropy energy barrier 
from 44 kT to 36 kT at room temperature. We assume 
that a permanent magnet will be employed to produce 
the bias field and thus will not affect the energy dissi- 
pated during switching. 
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D. Energy Dissipation 

The energy dissipated during switching has two com- 
ponents: (1) the energy dissipated in the switching cir- 
cuit that applies the stress on the nanomagnet by gener- 
ating a voltage, and (2) the energy dissipated internally 
in the nanomagnet because of Gilbert damping. We will 
term the first component 'CV 2 ' dissipation, where C and 
V denote the capacitance of the piezoelectric layer and 
the applied voltage, respectively. If the voltage is turned 
on or off abruptly, i.e. the ramp rate is infinite, then 
the energy dissipated during either turn on or turn off 
is (1/2)CV 2 . However, if the ramp rate is finite, then 
this energy is reduced and its exact value will depend 
on the ramp duration or ramp rate. We calculate it fol- 
lowing the same procedure described in Ref. Q. The 
second component, which is the internal energy dissipa- 
tion Ed, is given by the expression J Q r Pd(t)dt, where r 
is the switching delay and Pd(t) is the power dissipated 
during switchin g 21 i 22 



Distribution of e under thermal fluctuations 



Pd(t) 



a\j\ 



(l + a 2 )M v 



|T B (t) + T M (t)| 



(18) 



We sum up the power Pd(t) dissipated during the en- 
tire switching period to get the corresponding energy 
dissipation E d and add that to the 'CV 2 ' dissipation in 
the switching circuit to find the total dissipation E to tai- 
The average power dissipated during switching is simply 
E d /r. 

There is no net dissipation due to random thermal 
torque since the mean of the random thermal field is 
zero. However, that does not mean that the tempera- 
ture has no effect on either Ed or the 'CV 2 ' dissipation. 
It affects Ed since it raises the critical stress needed to 
switch with non-zero probability and it also affects the 
stress needed to switch with a given probability. Fur- 
thermore, it affects 'CV 2 ' because V must exceed the 
thermal noise voltag to prevent random switching due 
to noise. In other words, we must enforce CV 2 > kT. 
For the estimated capacitance of our structure (2.6 fF), 
this translates to V > 1.3 mV. 



III. SIMULATION RESULTS AND DISCUSSIONS 

In our simulations, we consider the magnetostrictive 
layer to be made of polycrystalline Terfenol-D that has 
the following material properties - Young's modulus 
(Y): 8xl0 10 Pa, magnetostrictive coefficient ((3/2)A s ): 
+90xl0 -5 , saturation magnetization (M s ): 8xl0 5 A/m, 
and Gilbert's damping constant (a): 0.1 (Rcfs. |24H27| ). 
For the piezoelectric layer, we use lead-zirconate-titanate 
(PZT), which has a dielectric constant of 1000. The PZT 
layer is assumed to be four times thicker than the magne- 
tostrictive layer so that any strain generated in it is trans- 
ferred almost completely to the magnetostrictive layer—. 
We will assume that the maximum strain that can be 
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FIG. 2. Distribution of polar angle 9i n ui a i and azimuthal an- 
gle (^initial due to thermal fluctuations at room temperature 
(300 K) when a magnetic field of flux density 40 mT is applied 
along the out-of-plane hard axis (+a;-direction). (a) Distribu- 
tion of polar angle 6 knitial at room temperature (300 K). The 
mean of the distribution is 173.7°, and the most likely value 
is 175.8°. (b) Distribution of the azimuthal angle 4>i n itiai due 
to thermal fluctuations at room temperature (300 K). There 
are two distributions with peaks centered at ~65° and ~295° . 



generated in the PZT layer is 500 ppm3£, which would 
require a voltage of 111 mV because d^i — l.^x 10~ 10 m/V 
for PZT—. The corresponding stress is the product of the 
generated strain (500 x 10 -6 ) and the Young's modulus 
of the magnetostrictive layer. Hence, 40 MPa is the max- 
imum magnitude of stress that can be generated on the 
nanomagnet. 

We assume that when stress is applied to initiate 
switching, the magnetization vector starts out from near 
the south pole (0 ~ 180°) with a certain (6 'initial, 4>initial) 
picked from the initial angle distributions at the given 
temperature. Stress is ramped up linearly and kept con- 
stant until the magnetization reaches the plane defined 
by the in-plane and the out-of-plane hard axis (i.e. the 
x-y plane, 9 = 90°). This plane is always reached sooner 
or later since the energy minimum of the stressed mag- 
net in 0-space is at 9 = 90°. Thermal fluctuations can 
introduce a spread in the time it takes to reach the x-y 
plane but cannot prevent the magnetization from reach- 
ing it ultimately if the stress is so large that the energy 
minimum at 9 — 90° is more than a few kT deep. 
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FIG. 3. Percentage of successful switching events among the 
simulated switching trajectories (or the switching probabil- 
ity) at room temperature in a Terfenol-D/PZT multiferroic 
nanomagnet when subjected to stress between 10 MPa and 
30 MPa. The stress at which switching becomes ~100% suc- 
cessful increases with ramp duration. For large ramp duration 
(120 ps) or slow ramp rate, ~100% switching probability is 
unachievable. 



As soon as the magnetization reaches the x-y plane, 
the stress is ramped down at the same rate at which 
it was ramped up, and reversed in magnitude to facili- 
tate switching. The magnetization dynamics ensures that 
8 continues to rotate towards 0° with very high proba- 
bility. When 8 becomes < 5°, switching is deemed to 
have completed. A moderately large number (10,000) 
of simulations, with their corresponding (6 initial, 4>initial) 
picked from the initial angle distributions, are performed 
for each value of stress and ramp duration to generate 
the simulation results in this paper. 

Fig. [2] shows the distributions of initial angles 6 knitial 
and (f>initiai in the presence of thermal fluctuations and an 
applied bias magnetic field along the +a;-direction. The 
latter has shifted the peak from the easy axis (8 = 180°). 
In Fig. 2(b) the <j)mitiai distribution has two peaks and 



resides mostly within the interval [-90°, +90°] since the 
bias magnetic field is applied in the +a;-direction. Be- 
cause the magnetization vector starts out from near the 
south pole (8 ~ 180°) when stress is turned on, the effec- 
tive torque on the magnetization — ItI / (l + v* 2 ) M x H 
due to the +a:-directed magnetic field is such that the 
magnetization prefers the ^-quadrant (0°,90°) slightly 
over the ^-quadrant (270°, 360°), which is the reason for 
the asymmetry in the two distributions of ^initial ■ Conse- 
quently, when the magnetization vector starts out from 
8 ~ 180°, the initial azimuthal angle ^initial 1S a lit- 
tle more likely to be in the quadrant (0°,90°) than the 
quadrant (270°, 360°). 

Fig. [3] shows the switching probability as a function 
of stress for different ramp durations (60 ps, 90 ps, 120 
ps) at room temperature (300 K). The minimum stress 
needed to switch the magnetization with ^100% prob- 
ability at K is ^5 MPa, but at 300 K, it increases to 
~14 MPa for 60 ps ramp duration and ~17 MPa for 90 ps 
ramp duration. At low stress levels, the switching prob- 
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FIG. 4. The thermal mean of the switching delay (at 300 K) 
versus stress (10-30 MPa) for different ramp durations (60 ps, 
90 ps, 120 ps). Switching may fail at low stress levels and also 
at high stress levels for long ramp durations. Failed attempts 
are excluded when computing the mean. 



ability increases with stress, regardless of the ramp rate. 
This happens because a higher stress mitigates the detri- 
mental effects of thermal fluctuations more when magne- 
tization reaches the x-y plane and thus conducive to more 
success rate of switching. This feature is independent of 
the ramp rate. 

Once the magnetization vector crosses the x-y plane 
(i.e. in the second half of switching), the ramp rate be- 
comes important. Now, the stress initially applied to 
cause switching becomes harmful and impedes switching. 
That happens because it causes the energy minimum to 
be located at 8 = 90° , which will make the magnetization 
backtrack towards this location during the second half. 
This is why stress must be removed or reversed imme- 
diately upon crossing the x-y plane so that the energy 
minimum quickly moves back to 8 = 0°,180°, and the 
magnetization vector rotates towards 8 = 0°. If the re- 
moval rate is fast, then the success probability remains 
high since the harmful stress does not stay active long 
enough to cause significant backtracking. However, if the 
ramp rate is too slow, then significant backtracking oc- 
curs whereupon the magnetization vector returns to the 
x-y plane and thermal torques can subsequently kick it 
to the starting position at 9 ~ 180°, causing switching 
failure. That is why the switching probability drops with 
decreasing ramp rate. 

The same effect also explains the non-monotonic stress 
dependence of the switching probability when the ramp 
rate is slow. During the first half of the switching, when 
8 is in the quadrant [180°, 90°], a higher stress is helpful 
since it provides a larger torque to move towards the x-y 
plane, but during the second half, when 8 is in the quad- 
rant [90°, 0°], a higher stress is harmful since it increases 
the chance of backtracking. These two counteracting ef- 
fects are the reason for the non-monotonic dependence of 
the success probability on stress in the case of the slowest 
ramp rate. 

Fig. 0] shows the thermally averaged switching delay 
versus stress for different ramp durations. Only success- 
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FIG. 5. The standard deviations in switching delay versus 
stress (10-30 MPa) for 60 ps ramp duration at 300 K. We 
consider only the successful switching events in determining 
the standard deviations. The standard deviations in switching 
delay for other ramp durations are of similar magnitudes and 
show similar trends. 



ful switching events are counted here since the switching 
delay will be infinity for an unsuccessful event. For a 
given stress, decreasing the ramp duration (or increas- 
ing the ramp rate) decreases the switching delay be- 
cause the stress reaches its maximum value quicker and 
hence switches the magnetization faster. For ramp du- 
rations of 60 ps and 90 ps, the switching delay decreases 
with increasing stress since the torque, which rotates the 
magnetization, increases when stress increases. How- 
ever, for 120 ps ramp duration, the dependence is non- 
monotonic, because of the same reasons that caused the 
non-monotonicity in Fig. 3. Too high a stress is harmful 
during the second half of the switching since it increases 
the chances of backtracking. Even if backtracking can be 
overcome and successful switching ultimately takes place, 
temporary backtracking still increases the switching de- 
lay. 

Fig. [5] shows the standard deviation in switching delay 
versus stress for 60 ps ramp duration. At higher values of 
stress, the torque due to stress dominates over the ran- 
dom thermal torque that causes the spread in the switch- 
ing delay. That makes the distribution more peaked as 
we increase the stress. 

Fig.[B]shows the thermal mean of the total energy dissi- 
pated to switch the magnetization as a function of stress 
for different ramp durations. The average power dissi- 
pation (Etotai/T~) increases with stress for a given ramp 
duration and decreases with increasing ramp duration for 
a given stress. More stress requires more 'CV 2 ' dissipa- 
tion and also more internal dissipation because it results 
in a higher torque. Slower switching decreases the power 
dissipation since it makes the switching more adiabatic. 
However, since the switching delay curves show the op- 
posite trend (see Fig.H]), the energy dissipation curves in 
Fig. |5] exhibit the cross-overs. 

Fig.[7]shows the 'CV 2 ' energy dissipation in the switch- 
ing circuitry versus stress. Increasing stress requires in- 
creasing the voltage V, which is why the 'CV 2 ' energy 



FIG. 6. Thermal mean of the total energy dissipation versus 
stress (10-30 MPa) for different ramp durations (60 ps, 90 ps, 
120 ps). Once again, failed switching attempts are excluded 
when computing the mean. 




FIG. 7. The l CV 2 ' energy dissipation in the external circuit 
as a function of stress for different ramp durations. 



dissipation increases rapidly with stress. This dissipation 
however is a small fraction of the total energy dissipa- 
tion (< 15%) since a very small voltage is required to 
switch the magnetization of a multiferroic nanomagnet 
with stress. The 'CV 2 ' dissipation decreases when the 
ramp duration increases because then the switching be- 
comes more 'adiabatic' and hence less dissipative. This 
component of the energy dissipation would have been 
several orders of magnitude higher had we switched the 
magnetization with an external magnetic field^S or spin- 
transfer torqueii. 

Fig. |8] shows the delay and energy distributions in the 
presence of room-temperature thermal fluctuations for 15 
MPa stress and 60 ps ramp duration. The high-delay tail 
in Fig. 8(a) is associated with those switching trajectories 
that start very close to 8 = 180° which is a stagnation 
point. In such trajectories, the starting torque is vanish- 
ingly small, which makes the switching sluggish at the 
beginning. During this time, switching also becomes sus- 
ceptible to backtracking because of thermal fluctuations, 
which increases the delay further. Since the energy dissi- 
pation is the product of the mean power dissipatio n and 
the switching delay, similar behavior is found in Fig. |8(bJ[ 

Fig. IHl shows two examples of switching dynamics when 
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FIG. 8. Delay and energy distributions for 15 MPa applied 
stress and 60 ps ramp duration at room temperature (300 
K). (a) Distribution of the switching delay. The mean and 
standard deviation of the distribution are 0.44 ns and 83 ps, 
respectively, (d) Distribution of energy dissipation. The mean 
and standard deviation of the distribution are 184 kT and 15.5 
kT at room temperature, respectively. 



FIG. 9. Temporal evolution of the polar angle 6(i) and az- 
imuthal angle (j>(t) for 10 MPa applied stress and 60 ps ramp 
duration. Simulations are carried out for room temperature 
(300 K). (a) Magnetization switches successfully, (b) Mag- 
netization fails to switch and backtracks towards the initial 
state. 



the applied st ress i s 10 MPa and the ramp duration is 
60 ps. In Fig. 9(a)| magnetization switches successfully. 
Thermal fluctuations cause the ripples because of tem- 
porary backtracking but 8 switches from ^180° to ^0° 
finally. Note that despite appearances, 4> is not changing 
discretely. When it crosses 360°, it re-enters the quad- 
rant [0°,90°], which is why it appears as if there is a 
discrete ju mp in the value of 4> in Fig. [9] On the other 
hand, Fig. |9(b)] shows a failed switching dynamic. Here, 
the magnetization backtracks towards 8 = 180° and set- 
tles close to that location, thus failing in its attempt to 
switch. This happened because of the coupled 6-(j> dy- 
namics that resulted in a misdirected torque when the 
magnetization reached the x-y plane. This kind of dy- 
namics has been explained in Ref. |T|| 



IV. CONCLUSIONS 

We have theoretically investigated stress-induced 
switching of multiferroic nanomagnets in the presence 
of thermal fluctuations. The room-temperature thermal 



average of the energy dissipation is as small as ^200 
kT while the thermal average of the switching delay is 
~0.5 ns with a standard deviation less than 0.1 ns. This 
makes strain-switched multiferroic nanomagnets very at- 
tractive platforms for implementing non- volatile memory 
and logic systems because they are minimally dissipative 
while being adequately fast. Our results also show that 
a certain critical stress is required to switch with ^100% 
probability in the presence of thermal noise. The value of 
this critical stress increases with decreasing ramp rate un- 
til the ramp rate becomes so slow that ^100% switching 
probability becomes unachievable. Thus, a faster ramp 
rate is beneficial. The energy dissipations and switching 
delays are roughly independent of ramp rate if switching 
is always performed with the critical stress. 
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